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Abstract 

The Doppler effects for the reflected wave from a moving media are systemically analyzed in this 
paper. The theoretical formula for the Doppler shift in the left-handed material, which is described 
by Drude's dispersion model, is presented. This formula is examined by first-principles numerical 
experiments, which are in agreement with the theoretical results. 
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I. INTRODUCTION 



The Doppler effect is a well-known phenomenon that the frequency of a wave emitting 
from a moving object is shifted from the source frequency, originally discovered by Christian 
Doppler in 1843 y], Conventional Doppler shift tells us that the measured frequency 
will increase when a source and observe approach each other. In virtue of this effect, an 
immense variety of novel applications are widely established including weather and aircraft 
radar system, blood flow measurement, and discovery of new planets and binary stars. The 
inverse Doppler effect refers to frequency shifts that are in the opposite sense to those 
mentioned above. It has been observed in a dielectric-loaded stripline cyclotron-resonance 
maser oscillator [_3jj, the moving shock wave propagating through a photonic crystals ^| 
and the moving boundary in an electrical transmission line Q]. In fact, back to 1968, the 
inverse Doppler effect has been already predicted to occur in a left-handed material (LHM) 
with simultaneously negative permittivity and permeability 6:]. Pendry et al. discovered a 
microwave-plasma thin-wire structure exhibiting negative permittivity below the electronic 
plasma frequency [?J and a magnet-free split-ring resonator structure exhibiting negative 



permeability below the magnetic plasma frequency 



. With the suggestion of Pendry's 



theory, the first experimental verification of a LHM was demonstrated |9j, which has recently 



received much attention in the literature 
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13j,ll4j. However, up to our knowledge, 



inverse Doppler effect in a left-handed material has not been demonstrated experimentally 
yet, not even by numerical experiments. 

In this paper, we will study Doppler effects for the reflected wave from a moving media and 
present theoretical formula for the Doppler shift, which will be examined by first-principles 
numerical experiments using the finite difference time domain (FDTD) method. 



II. MODEL AND THEORETICAL FORMULA 



Consider a plane electromagnetic wave propagation in medium / and incident at an angle 
9\ on an idealized moving medium II. The incident wave have an angular frequency uj\ 
and a propagation vector kx (In a LHM, the direction of group velocity, same as that of 
energy flux, is often antiparallel to phase velocity). To calculate Doppler shift effectively for 
the reflected wave in such an inertial system, a primed coordinate frame embedded in the 
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moving medium is introduced. The primed frame moves with the same velocity v as the 
moving medium relative to the unprimed (original) coordinate system. Here we keep the 
axes of the primed coordinate system parallel to corresponding axis of the unprimed system. 
Therefore, we can treat the space components of both the primed and unprimed system as 
vectors in a common three-dimensional space. 

If the primed frame is moving with a velocity v with respect to the unprimed, one can 
relate r' and t' with r and t through a pure Lorentz transformation (no rotation of the space 

n 

axes) expressed in three-dimensional notation (see Podolsky and Kunz |15l|). 

r = &(v)r - yvtj = y(t - (l/c 2 )v T r), (1) 

where 

7 = (1 - v 2 /c 2 Y l '\ *(v) = 1 + (7 - l)vv T , (2) 

and where v = v/|v| is the unit vector in the direction of v and v T is the transpose of v. 
Consider a 4- vector expressed in the Minkowski notation by (k, iuj/c). Based on the theorem 
of invariant phase of an electromagnetic wave for the transformation, one can easily obtain 

k' = $(v)k - 7 v(cu/c 2 ), J = y(uj - v T k), (3) 

which is the Doppler shift in going to the primed frame. 

To obtain a solution to the problem posed above it is convenient to firstly consider in the 
primed frame, in which the medium II is at rest while the background moves with a velocity 
— v. In the primed frame there is no Doppler shifts at the interface, we can set 

uj[ = J 2 (4) 

as the required generalization, where uj[,u) 2 is the frequency of incident and reflected wave 
in the primed frame. Using the inverse of Eq. (jSJ), one can replace Eq. (j3J) and then obtain 

vigi = u 2 g 2 , (5) 

where 

9i = 1 - v T ki/V (6) 
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Finally, one can obtain the reflected frequency in respects of the incident frequency as 

Ul 2 1 + Ti\V cos 9i/c 



(7) 



ui 1 — n'xV cos 6-i/c 

where 9\ and 62 denotes the incident and reflected angle, and where n\ and n[ denotes the 
refractive index of the background for the incident and reflected wave, respectively. In the 
case of moving media, the law of reflection and Snell's law of refraction is no longer suitable 
at the interface. Moreover, since the LHM is often a dispersive material, the refractive index 
(rii) of the background for the reflected wave is not the same as that of the background for 
the incident wave (ni) due to the Doppler shift. 

As a LHM is always a dispersive media, it can be described by the following Debye model, 



e, = eoo 1 + - f— — , (8) 

Vt = ^oo[l + — U 2 \ . ) ■ (9) 

When = 1, a>& = 0, and ui a = u p , it will be simplified to be the following Drude's lossless 
dispersive model for the relative permittivity and permeability of the LHM 

uj 2 



e r (u) = 1~J, (10) 

A*rH = l-4« ( H ) 

where u p is the plasma frequency or the magnetic resonance frequency. The permittivity and 
permeability take negative values for frequencies below u p , which satisfy e r (ui) = ft r (^i) = 
— 1 at frequency u\ = tu p /\/2. Combining Eq. ((7J) with the definition of negative refractive 



index in the LHM (n = —Je r (u)[i r (u)), we can obtain the Doppler shift theoretically for 
the reflected wave 

^i(l + mv/c) + J[uji{1 + niv/c)} 2 - 4^(1 _ v /c)v/c 

" 2 = W^Jc) ' (12) 

with m = 1 — ujpi/uf. Here n\ and u p i is the refractive index of the background for 
the incident wave and the plasma or magnetic resonance frequency for the background, 
respectively. Using the above formula, the Doppler shift can be easily obtained by an 
analytical way if we know the incident frequency and the parameters in each medium. 
Absorption will not exit for this idealized model for the LHM. Moreover, in this paper, one 
is interested with the frequency shift of the reflected wave, not the intensity of the reflected 
waves. 
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FIG. 1: Doppler shifts for the reflected wave from the moving medium in air, which is a conventional 
material with a refractive index of 11,2 = 3. The open circle markers represent the values predicted by 
the Doppler formula, and the asterisk markers represent the numerically simulated values obtained 
by the FDTD method. 

III. NUMERICAL SIMULATIONS AND RESULTS 



Next, we will give some numerical examples to examine the Doppler formula 



by first- 



ly. For 



principles numerical experiments using the finite difference time domain method 
simplicity, we consider the .E-polarized electromagnetic wave incident normally on the mov- 
ing media, which moves toward the source. It is similar for the case of the media moving 
away from the source. We perform firstly a FDTD simulation to study the case of two 
conventional materials. The background is air and the refractive index of the right moving 
medium is n 2 = 3. A line source of continuous wave along y direction is placed at the left 
side of the moving media, as well as the detector. Since we only consider the electromagnetic 
wave incident normally on the interface, it is naturally chosen the periodic condition in the y 
direction and the perfect matched layers (PMLs) [18] in the x direction as numerical bound- 



ary treatments. Due to stability of the FDTD algorithm (At < l/(cy 1/ Ax 2 + 1/ Ay 2 )), the 
moving interface always moves less than Ax for each time step. In all our simulations, a 
technique of parameter approximation based on the first order Lagrange approximation is 
used. The theoretical Doppler formula for the conventional material is the same as that for 
the case of the LHM, except frequency-independent refractive index for conventional mate- 
rials in Eq. (jZJ . Figure d shows the Doppler shifts for the reflected wave from the moving 
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media. In Fig. ^ the open circle markers represent the values predicted by the Doppler 
formula, and the asterisk markers represent the numerically simulated values obtained by 
FDTD method. It can be seen from Fig. [J that the results of our simulations are in good 
agreement with the theoretical results. From Fig. Q one also finds that the frequency of 
reflected wave becomes larger for the media moving towards the source, which is referred as 
normal Doppler shift (blueshifted). 
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FIG. 2: Doppler shifts for the reflected wave from the LHM moving in air. The refractive index 
of the LHM for the incident frequency is —1, matched with the background (air). The open circle 
and asterisk markers have the same meaning as those in FIG. 1. 



For the LHM described by Drude's model, as a special case, we firstly consider the system 
with a LHM moving in air. The FDTD method for the Drude's model is discussed in detail 
in Q. We choose . = ^ whet, ^ is the plastna fluency ot tnagnetic _ 
frequency of the LHM. At that time, the refractive index of the LHM for the incident 
frequency (ui) satisfies n 2 = — 1, matched with air. It has been shown that light can go 
through such an stationary air-LHM interface without reflection for any incident angle [20]. 
However, this is only correct for the case of a stationary interface since the field boundary 
conditions at moving interface are quite different with those at the stationary case . The 
Doppler shift for the reflected wave is governed by Eq. (J7J) for the electromagnetic wave 
incident normally on the moving interface. Results for the Doppler shifts for the reflected 
wave are shown in Fig. |2j The values predicted by the Doppler formula (Eq. (|12j) ) are 
presented by open circle markers, and simulation results obtained by the FDTD method 
are shown by asterisk markers. From Fig. El one finds that simulation results are in good 
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FIG. 3: The effective refractive index of the background for the reflected wave as the velocity of 
the moving medium increases. 

agreement with the theoretical values, which are also consistent with the results in Fig. 
It can be understood from Eq. (JZj) that the Doppler shift for the reflected wave only relates 
with velocity of the moving media and electromagnetic property of the background, no any 
relation with the moving media. 

As stated in the introduction, the reverse Doppler effect arises due to the negative index 
of refraction. The definition of Doppler shift tells us that when a source moves inside the 
LHM at the initial stationary frame, the abnormal Doppler effect will be observed. Here we 
consider the system with a LHM moving in the background with another LHM. We choose 
u)\ = L0 p / y/2, ujpi = L0 P , and uj p2 = V^^p, which means that the refractive index of the 
background and the moving media for the incident frequency is —1 and —3, respectively. 
Due to Doppler shift for the reflected wave, the refractive index of the background for the 
reflected wave will change. Figure El shows that the refractive index of the background 
for the reflected wave varies from —1.59 to —3.60 (Note that the refractive index of the 
background for the incident frequency is —1). The Doppler shift for the reflected wave from 
the moving media are shown in Fig. HJ in which the circle and asterisk markers represent 
the theoretical values (obtained by Eq. ([12ft ) and simulation results (obtained by the FDTD 
method), respectively. Although the simulation results have a difference with theoretical 
values, the numerical errors are always less than 1.5%. The errors may be caused by the 
numerical calculation of the FDTD method. It can be seen from Fig. EJthat, compared with 
the results in Figs. Hand El, the Doppler shift for the reflected wave is actually an inverse 
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FIG. 4: Doppler shifts for the reflected wave from a LHM moving in a background with another 
LHM. The refractive index of the background and the moving media for the incident wave is — 1 
and —3, respectively. The open circle and asterisk markers have the same meaning as those in 
FIG. 1. 

Doppler effect (redshifted) . It is also in agreement with the prediction in the LHM. 



IV. CONCLUSION 



We have studied systemically the Doppler effect for the reflected wave from a moving 
media and obtained theoretical formula for the Doppler shift in the LHM, which is described 
by Drude's dispersive model. We have performed first-principles numerical experiments to 
examine the Doppler shifts. It has been shown that the results obtained by our theoretical 
formula are in good agreement with those obtained by numerical experiments. Inverse 
Doppler effect is confirmed in the left-handed material. 

This work was supported by the Swedish Foundation for Strategic Research (SSF) on IN- 
GVAR program, the SSF Strategic Research Center in Photonics, and the Swedish Research 
Council (VR) under project 2003-5501. 
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